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ABSTRACT 

One of the factors that influences the fidelity of an aircraft digital simulation is the 
sampling rate. As the sampling rate is increased, the calculated response of the discrete 
representation tends to coincide with the response of the corresponding continuous system. 
Because of computer limitations, however, the sampling rate cannot be increased indefinitely. 
Moreover, real-time simulation requirements demand that a finite sampling rate be adopted. 

In view of these restrictions, a 9tudy was undertaken to determine the influence of sampling 
rate on the response characteristics of a simulated aircraft describing short-period oscil- 
lations. Changes in the calculated response characteristics of the simulated aircraft 
degrade the fidelity of the simulation. In the present context, fidelity degradation is 
defined as the percentage change in those characteristics that have the greatest influence 
on pilot opinion: short-period frequency w, short-period damping ratio £, and the product 
Cu>. To determine the influence of the sampling period on these characteristics, the equa- 
tions describing the response of a DC-8 aircraft to elevator-control inputs were used. The 
results indicate that if the sampling period is too large, the fidelity of the simulation 
can be degraded. For example, it was determined that the motion characteristics of the simu- 
lated aircraft, represented by t and £u>, as well as the corresponding fidelity, decreased 
linearly with sampling period and vanished when the period reached 640 ms. 

INTRODUCTION 

The response of a dynamical system to control inputs may be stable or unstable, depending on 
the characteristics of the system. When the response is being calculated by a digital com- 
puter, the dynamical characteristics are modified by the need to use a discrete representa- 
tion of the system. This involves the choice of a sampling rate, which influences the 
fidelity of the simulation and determines the feasibility of simulating large aeronautical 
systems, such as helicopters. 

One way to determine the influence of a specific parameter on the response characteristics 
of a continuous dynamical system is to map the eigenvalues on the complex plane and note 
their migrations, in each mode, a9 the parameter is varied (1). When the loop is closed 
around an open-loop system, and the feedback gain is varied, the conventional root loci are 
obtained (2). These loci reveal the stability characteristics of the system, as the fre- 
quency and the damping ratio vary with changes in the feedback gain. 

By using the same procedure for a discrete representation of the system, the influence of 
sampling rate on the response characteristics can be determined. Whereas a linear, continu- 
ous dynamical system is described in terras of the Laplace transform, the discrete represen- 
tation requires a formulation in terms of the Z-transform. When the Z-transfer function 
corresponding to a given integration algorithm has been obtained as a function of the sam- 
pling period, and the discrete eigenvalues determined, these can be plotted on the complex 
z-plane to ascertain their location relative to the unit circle and thus determine the pos- 
sibility of instability. To facilitate the interpretation of results, the discrete eigen- 
values were transformed from the z-plane to the s-plane. By proceeding in this manner, and 
using the equations describing the response of a DC-8 aircraft to elevator-control inputs, 
the influence of sampling rate was determined. 
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ANALYSIS 


Equations of Motion 


The differential equation describing an aircraft's pitching motion, which is assumed to be 
linear and time invariant, is (from Ref. 3) 

/ d 4 Y , . d 3 Y . . d 2 Y , dY , , v \ /d 2 U . n dU , D 

W 3 ^ 2 ^ 1 * °7 ‘ W 1 * °7 

where Y is the pitch angle, U is the elevator-control input, and A^ , Bj are constants. 

The formulation is simplified by specifying a state variable X which satisfies the differ- 
ential equation 


d**X 


d J X 


d z X 


dX 


__ + A 3 _ + A 2 — + Al ^ + A 0 X - U 

The remaining state variables can then be assigned as follows: 

X - x 1 
*1 - x 2 

X 2 - X, 

x 3 - x„ 

X„ - -A 0 Xj - AjXj - a 2 x 3 - a 3 x„ + u 

Using matrix notation, these equations assume the conventional form; that is, 

X ■= AX + BU 


where 


and 


From Equation (4) , 


Y - C A X 

X - (X x X 2 X 3 X„) T 

0 10 

0 0 1 

0 0 0 

B T - (0 0 0 K) 

C T - (B 0 B x 1 0) 


X(t) - e A(t ’ to) X(t 0 ) + J e A(t " T) BU(T)dT 


( 2 ) 


(3) 


(A) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


When the input is sampled with a zero-order hold after the sampler, the input is held con- 
stant for one sampling period (see Figure 1) . 

During this period the limits of integration are t - (n + 1)T and t 0 * nT, and Equa- 
tion (4) becomes (from Ref. 4) 

X(n + IT) - e AT X (nT) + ( } e A(n+lT " T) dx^BU(nT) 
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Therefore, 


X(n + IT) - e AT X(nT) + A -1 (e AT - I)BU(nT) (11) 

By omitting the T from the arguments of X and U, this equation assumes a more compact 
form, namely, 

X(n + 1) • 4>x(n) + B^Cn) (12) 

where - e AT , and B D is the discrete form of the control distribution matrix. That is, 

B d - A -1 (♦ - I)B (13) 


It should be noted that Equation (13) is valid only if A is nonsingular. 

Z-Transforms and Z-Transfer Functions 

Assuming zero initial conditions, the Z-transform of the state vector is obtained from 
Equation (12) : 

X(Z) - (ZI - »)- l B D U(Z) (14) 

The corresponding Z-transform of the output is 

Y (Z) - C T (ZI - t) _1 BjjU(Z) (15) 

The Euler Integration algorithm usea only the first two terms of the transition matrix e A *: 

e AI ■ * f (I + AT) (16) 

Substituting this value into Equation (13) gives the Euler form of the control distribution 
matrix B D , 

B D - A -1 (I + AT - I)B - TB (17) 

Substituting this result into Equation (15) gives the Z-transform of the output correspond- 
ing to the Euler integration algorithm: 

Y (Z) - TC T [(Z - 1)1 - AT]” 1 BU(Z) (18) 

The Z-transfer function relating the output Y(Z) to the input U(Z) is 

Y(Z)/U(Z) - TC T [(Z - 1)1 - AT] _l B (19) 


By introducing matrix A from Equation (7) and the column and row vectors from Equations (8) 
and (9), respectively, the required form of the transfer function is obtained: 


where 


Y(Z) 

U(Z) 


xr 


/ z 2 + P X Z + p 0 \ 

\z k + q 3 z 3 + q 2 z 2 + q x z + q 0 / 


( 20 ) 


T - sampling period 
p 0 ■= B„T 2 - BJ + 1 
Pl - B lT - 2 

q 0 =■ AjI 1 * - AjT 3 + AjT 2 - AjT + 1 

q 1 - A 1 T 3 - 2A 2 T 2 + 3A 3 T - A 
q 2 - A 2 T 2 - 3A 3 T •»* 6 
q 3 “ a 3 t - 4 


RESULTS 


The following values of the elements of the system dynamics matrix A and the control dis- 
tribution matrix B were obtained from Teper (3). 
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A 0 - 0.07006953 


A x - 0.09712526 


A 2 - 2.6813872 K - 1.338 


A 3 - 1.700522 


B x - 0.5955 


B 0 - 0.032675 


When these values were substituted into Equation (20), the open-loop eigenvalues correspond- 
ing to the phugoid and the short-period modes were obtained as functions of the sampling 
period T. To determine the influence of the sampling rate on the stability of the short- 
period motion, the discrete eigenvalues were plotted as a function of the sampling period, 
and their locations relative to the unit circle were noted. The results are plotted in 
Figure 2. It is seen that as the sampling period is increased, the eigenvalues remain 
inside the unit circle until T - 640 ms. At this point they lie on the unit circle, and 
further increases in T cause them to move into the region of instability. 

In order to determine the changes in frequency and damping as functions of the sampling 
period, the discrete eigenvalues were transformed from the Z-plane to the s-plane. The 
migrations of the eigenvalues in the s-plane give a good indication of the influence of the 
sampling period on the frequency and the damping of the simulated system. The required 
transformation can be effected by using the equation of definition: 


Equation (21) was used to calculate the T-locus of open-loop eigenvalues for the following 
range of sampling periods: 


The results are plotted in Figure 3, where it is seen that the locus reaches the imaginary 
axis when the sampling period is 640 ms. This figure should be compared with Figure 2, where 
the locus of discrete eigenvalues reaches the unit circle for the same value of T. More- 
over, the data used to plot Figure 3 can be used to calculate the modal frequency and the 
damping ratio as functions of the sampling period. The T-locus of frequency is plotted in 
Figure 4, and the locus of the corresponding damping ratios is shown in Figure 5. When T 
is increased from 20 to 640 ms, it appears that the damping ratio decreases almost linearly 
from its initial value to zero. The undamped natural frequency is shown plotted as a func- 
tion of the damping ratio in Figure 6. The range of this curve indicates that as the sam- 
pling period is increased, unsatisfactory regions of the flying-qualities plot will be 
encountered (5). When simulation data become available, boundaries of acceptable flying 
qualities will be established. These boundaries assume the form shown in Figure 7, where 
acceptable and unacceptable flying qualities are related to the short-period frequency and 
to the short-period damping ratio (6). 

Because the response of a simulated dynamical system to control inputs is calculated by a 
digital computer, the motion characteristics of the discrete representation are functions of 
the sampling rate; consequently, simulation fidelity is degraded. In the present context, 
fidelity degradation is defined as the percentage change in those characteristics that have 
the greatest influence on pilot opinion: short-period frequency w, short-period damping 

ratio c» and the product For example, the fidelity loss owing to changes in the damp- 

ing ratio is defined as 


where f^ is the fidelity loss, £ c is the continuous damping ratio, and Gd the damping 
ratio of the discrete representation. The f^ term is plotted as a function of the sampling 
period in Figure 8. As might be expected, f^ has a functional form similar to that of c, 
as shown in Figure 5. 

The loss of fidelity that results from natural frequency changes is defined similarly: 


s - (log e Z)/T 


( 21 ) 


0.02 < T < 0.80 s 


( 22 ) 



(23) 



( 24 ) 
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where f w is the fidelity loss, oj c is the natural frequency of the continuous system, and 
is the natural frequency of the discrete representation. As shown in Figure 8, the 
influence of natural frequency changes on system fidelity, as measured by Equation (24), is 
not so pronounced as the influence of the damping ratio. The third measure considered 
assesses the loss of fidelity owing to the characteristic Ctu. The loss of fidelity caused 
by changes in £u> is defined as 


f 


Co) 


C c w c ~ Vb l 

^c w c 


(25) 


where again the subscripts c and D denote continuous and discrete representations, respec- 
tively. This characteristic brings out the influence of the real component of the natural 
frequency vector on simulation fidelity. Since that component determines the damping ratio, 
the loss of fidelity, as defined by Equation (25), is seen to follow the same trend as f^. 
The results are plotted in Figure 9. 


CONCLUSIONS 

One of the factors that influences the fidelity of an aircraft simulation is the sampling 
rate. As the sampling rate is increased, the calculated response to control inputs, of the 
discrete representation, tends to coincide with the response of the continuous system. 

Because of computer limitations, however, the sampling rate cannot be increased indefinitely, 
and simulation fidelity is degraded. In the present context, fidelity degradation is 
defined as the percentage change in those characteristics that have the greatest influence 
on pilot opinion: short-period frequency w, short-period damping ratio 4, and the product 

too. In order to determine the influence of the sampling period on these characteristics, the 
equations describing the response of a DC-8 aircraft to elevator-control inputs were used. 
Results indicate that if the sampling period is too large, the fidelity of the simulation can 
be degraded significantly. For example, it was determined that the motion characteristics of 
the simulated aircraft, represented by £ and 4w, as well as the corresponding fidelity, 
decreased linearly with sampling period and vanished when the sampling period reached 640 ms. 
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Figure 1 — State Variable Sampled-Data System. 



Figure 2 — T-Locus for Solution of Equation (20) by Euler Integration. 
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Figure 3 — T-Locus of Open-Loop Poles for the Short-Period Mode. 
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Figure 4 — T-Locus of Natural Frequencies for the Short-Period Mode. 



Figure 5 — T-Locus of Open-Loop Damping Ratios for the Short-Period Mode. 
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Figure 6 — T-Locus of Undamped Natural Frequency versus Damping Ratio. 
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Figure 7 — Flying-Qualities Boundaries. 
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Figure 8 — Loci of Fidelity Measures for DC-8 During Landing Approach. 



Figure 9 — Locus of Fidelity Measures for DC-8 Aircraft in Approach Condition. 
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